Elliptic flow in high multiplicity proton-proton collisions at y/s =14 TeV as a signature 
of deconfinement and quantum energy density fluctuations. 
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At LHC extreme values of energy density will be reached even for proton-proton collisions. Such 
values of energy density may be large enough to generate a collective motion in the products of 
the collision, therefore generating effects such as elliptic flow. Using ideal 3+ ID hydrodynamical 
simulations, we show that elliptic flow can occur at least for top multiplicities p-p events at LHC 
and that the intensity of such effect is strongly related to quantum fluctuations in the initial proton 
energy distribution. 
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Introduction 

The application of hydrodynamical models to de- 
scribe hadronic collisions have a long history, dating 
almost from the beginning of the studies of hadronic 
interactions p], 0,0]. Since late '80s these approaches have 
been revived to describe the collective features of the dy- 
namics in heavy ion collisions at the big accelerators US 
[H and investigate the properties of quark-gluon plasma 
(QGP). The main idea behind such approach is that from 
the analysis of collective flow in the dynamics, if any, we 
may infer the properties of the allegedly deconfined par- 
tonic matter generated in heavy ion collisions. For this, 
the thermalization of the partonic matter should occur 
in a relatively short time and space scales. Until today 
this approach was limited to heavy ions collisions be- 
cause only in these collisions a large space-time domain, 
in which the energy density is sufficiently high to lead 
to the deconfinement and thermalization, is expected to 
emerge. Proton-proton (p-p) data were used as reference 
to analyze these heavy ion data. 

In p-p processes the system is considered too small and, 
so far, it was not expected to get any deconfinement nor 
thermalization domains. Thus, any differences in behav- 
ior between the p-p and heavy ions case may (should) be 
attributed to the dense medium effects. 

As mentioned above, one of the most important sig- 
natures of the creation of such a medium is the mani- 
festation of collective motions (flows) in the final state 
particles. The presence of any collective flow is a typi- 
cal behaviour of a fluid matter. There are many ways to 
identify collective flows. They appear as an anisotropy 
in angular distribution of particles and are referred to 
as anisotropic flow. In experimental studies, anisotropic 
flows are usually analyzed in terms of the Fourier expan- 



sion of the azimuthal particle distribution 

1 dN =1 + 2^1 cos(^) + 2v 2 cos (20) 
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where usually v\ is referred as direct flow and v 2 as el- 
liptic flow. 

It should be noted that these anisotropic parameters are 
not necessarily null even for free-streaming particles, but 
their dependence on quantities such as, for example, im- 
pact parameter and kinematic variables will reflect the 
collective features of the dynamics. 

LHC accelerator, that is starting in Geneva, will reach 
extreme values of energy both for heavy ions and pro- 
tons, the latter being accelerated up to y/s = 14 TeV. 
It is possible that under these extreme conditions, also 
in proton-proton collisions the energy density will reach 
values high enough to allow partonic matter to survive 
in a deconfined state for a time long enough to thermal- 
ize and show some collective behaviour. In particular, 
when we consider quantum fluctuations in the incident 
states, such a situation may well occur and thus mani- 
fest in large multiplicity events. In this sense, collective 
features, if any, of the high multiplicity events may reveal 
interesting informations on the initial condition for high 
energy p-p collisions. 

The goal of this work is to estimate the elliptic flow for 
top multiplicity proton-proton events in the framework 
of ideal hydrodynamics and investigate if it will be pos- 
sible to perform such analysis at the LHC. If viscosity is 
present, it has effects on collective flow and to determine 
the viscosity of the QGP fluid is presently one of the 
important problems. However, these questions are still 
under investigation, and in the present work, we only 
consider the case of an ideal fluid to avoid additional 
ambiguities associated to viscous theories [7]. 
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I. EOS AND INITIAL CONDITIONS 

For simplicity, we use a factorized initial energy density 
profile into its longitudinal and transverse parts, as pro- 
posed by [8]. For the transverse part, we use an energy 
density profile as a simple superposition of those of col- 
liding particles as is usually done for heavy ions collisions 
[9]. Thus we use the formula 



e (x,y,b) = 
k(y/Z)T A (x + l,y) 



k(yfi)T B (x--,y) 
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where A and B are the number of nucleons, a is the 
nucleon-nucleon cross section and k a collision energy de- 
pendent constant. In the case of proton-proton collisions 
all this is reduced to the simple expression 

e t=0 (z, y, b) = 2kaT A (x + ~,y)T B (x -^y). (3) 

For the proton energy density T^m, we use a Gaus- 
sian distribution with a top value eo and width 
a p = 0.875 fm. 



We have chosen entropy values between 5400- 
5500 fm -3 in order to reproduce the top multiplicities ex- 
pected at LHC (HIJING simulations) for such collisions. 
It is very important to focus on very high multiplicity 
events for two reasons. The first is that proton-proton 
collisions on average have very low multiplicity. At LHC 
energies for proton-proton collision we expect dN / dy w 6 
while the top multiplicity can be larger then one thou- 
sand particles. To extract any information on "collec- 
tivity", the events producing a large enough amount of 
particles are interesting. The second reasons is that in 
low multiplicity events also non-flow correlations may be- 
come strong to override the elliptic flow, even if it exists. 
Fixing the total multiplicity, we further determine the 
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T A (x,y) = e e~ 



(4) 



As the convolution of two Gaussian is still a Gaussian, 
our initial transverse energy density is 



e t =o(x,y,b) = Ke 



(5) 



where all the constants (fc,cr...) are compacted together 
in the new constant K. It must be noticed that the result 
of eq. ([5]) is completely symmetric in the transverse plane 
for any value of b. We will discuss this point in more 
detail later in section ITU 

For the longitudinal dimension, we proceed as in [8] 
using the equation: 



eiong(Vs) = e *%0(y\ 



'be am 



Vs)- 



(6) 



We are left with two free parameters, the constant K 
and the longitudinal width r] g . We can express the total 
entropy as a function of these parameters, and so relate 
them to the particle multiplicity, as shown in Fig. (pQ). 
Remarkably, the results on the integrated value of V2 are 
almost not affected by the choice of these parameters, 
provided that they keep the total entropy constant. Fur- 
thermore, it depends explicitly on the shape of the impact 
region, as is discussed later. Nonetheless, the momentum 
dependence of ^(pt) is strongly affected by those param- 
eters (cf. sec. HI]) . 



FIG. 1: (Color online) entropy dependence on the parameters 
rjg and K. The solid line shows the entropy region reproducing 
the top LHC multiplicities. 

parameter rj g by the expected mid-rapidity multiplicity 
of particles. This choice leads to the values rj g = 1.3 
and K = 1150 GeV. In order to test the momentum 
dependence of V2 on the parameter choice, other two 
sets of parameters are used in the paper: the first one 
(K = 370 GeV and rj g = 3.1) with very low top en- 
ergy density but extremely broad distribution in the rjr 
dimension, and the second one (K = 3680 GeV and 
rjg = 0.6) with a very high initial energy density but 
almost confined in the transverse plane. 

In ideal hydrodynamics the elliptic flow is proportional 
to the initial eccentricity of the impact region [1Q| . 
In our case, there is a big difference from heavy ion col- 
lisions. Due to the Gaussian shape of the proton energy 
densities, the total energy density profile eq.([3]) is sym- 
metric, so that at the classical level we would not expect 
elliptic flow at all. It has already been shown that quan- 
tum effects on energy distribution are not much effective 
on large systems, such as Au-Au collisions [11] but their 
importance grows as the size of the system is reduced, like 
for Cu-Cu collisions [l2[. So for the very small systems 
expected in p-p collisions, quantum fluctuations should 
have a crucial effect on the event-by-event energy den- 
sity distribution of the created matter. We try to mimic 
such quantum fluctuations by simply cutting the shape 
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of the interaction region with an elliptic shape parame- 
terized by eccentricity e. We thus suppose that the initial 
energy density for the hydrodynamical evolution as 

e(f, v, e) = Ke~^ e~^0(y beam - tj s )6 1 - ^ - f, , 
with 

a = \/l ( 1-e )' 

& = \/f(l + e). 

According to our hypothesis, the parameter e measures 
the effect of quantum fluctuation in the formation of the 
initial condition for ultra-high energy p-p collisions. 

As for the equation of state (EoS), we use the one 
shown in Fig.© (solid line). This EoS is obtained by 
smoothly connecting the equations of state of the lattice 
QCD result in [l3j] and an hadron resonance gas near 
T c = 200 MeV. The phase transition is a cross-over. To 
see the effect of the EoS on our results, we also test the 
EoS with a steeper connection between these two phases 
(dashed). After the initial entropy distribution is gen- 
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FIG. 2: (Color online) Behaviour of squared speed of sound 
and pressure over T 4 as functions of temperature for the two 
different equations of state. 

erated, the time evolution of the system is followed us- 
ing the Smoothed Particle Hydrodynamic (SPH) formal- 
ism in hyperbolic coordinate system To calculate 
the elliptic flow, we perform the commonly used sudden 
freeze-out through the Cooper-Frye procedure: when a 
fluid element crosses the surface defined by the tempera- 
ture T = Tf all interactions are assumed small enough so 
that all particles can be consider as free. The freeze-out 
temperature is settled at the value Tf = 130 MeV. Tests 
with Tf = 120 MeV shows that our results are almost 
insensitive to the changes in the value of Tf. 

A full 3+ ID ideal fluid dynamics using our SPH code 
for p-p requires at least 30 000 total number of SPH 
particles with the kernel width h = 0.3 fm [l4|. With 
this condition, one collision event consumes ~8 hours of 
calculation with a desktop PC. The overall precision of 
calculation is monitored by the conservation of the total 
energy and momentum which are kept within the relative 
error of 10 -2 during the whole time evolution. 



II. RESULTS AND CONCLUSIONS 
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FIG. 3: Azimuthal thermal pions distribution for 100 high 
multiplicity events. Fit leads to V2 ~ 3%. 

In Fig. (|3]) we present the azimuthal distribution of 
thermal pions in the mid rapidity region y G [—1, 1] for 
an initial eccentricity of 17%. The results are consis- 
tent with an elliptic flow of 2.9%, so this indicate that 
an elliptic flow effect could be strong enough to be ex- 
perimentally measurable at LHC. In Fig. ([4]) we show 
the effect of the choice of the parameters on V2(pt) for 
different values of initial eccentricity. As expected, the 
elliptic flow effect gets stronger when most of the en- 
ergy density is available in the transverse plane and the 
integrated value of V2 is proportional to the initial ec- 
centricity [10[. We also verify that the behavior of V2 
integrated over the whole rapidity domain is almost the 
same as that of the mid-rapidity region. If the eccentric- 
ity e has a quantum origin in p-p collisions, we expect 
big fluctuations for the average value of V2 for different 
multiplicity bins in an experiment. Also, we should ex- 
pect big fluctuations even in the same bin for different 
events, as the origin of the eccentricity is not geometrical 
but stochastic. We also expect the elliptic flow value to 
be always greater than 0, as any fluctuation in any col- 
lision contributes with a positive value of V2- However, 
the collective flow can only become important for large 
multiplicity events. From our study, we get the elliptic 
flow of V2 ~ 3% if the induced eccentricity has average 
value of about 17%. We expect that, measurements of 
elliptic flow in top multiplicity events in p-p may give 
an important clue for the quantum fluctuation in energy 
density distribution at each event. 

The ALICE collaboration is developing several meth- 
ods to estimate non-flow contribution to V2 in p-p 
collisions [l5|. The most promising one seems to be the 
77-gap method, that consists in using the particles with 
pseudorapidity 77 > 3 to evaluate the reaction plane and 
those at I77I < 1 to evaluate V2. Other methods are, for 
example, the standard event-plane method [16] or defin- 
ing the reaction plane in the direction of the leading par- 
ticle at high pseudorapidity. Once non-flow correlations 
are removed, V2 can be measured. Unfortunately, there 
aren't yet any available estimations of the precision such 
methods can reach in removing no n- flow correlations, but 
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it is plausible that an effect of a few percent would be vis- 
ible in ALICE. 

As mentioned in the introduction, one thing that can af- 
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feet the V2 value is the presence of viscosity in the collec- 
tive motion, that has been neglected in our calculations. 
High value of viscosity can strongly affect the value of 
the elliptic flow produced in proton-proton collisions, as 
suggested in [17[. 
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FIG. 5: (Color online) V2(pt) for the two different EoS pre- 
sented in fig. ( [2) for eccentricity 17%, K=1150 MeV and 
rj g 1.3. 
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FIG. 4: (Color online) V2(pt) for an initial eccentricity of 17% 
(up), 8.6% (center) and 26% (bottom) for the three different 
set of parameters 
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